The high temperature expansion of the classical XYZ chain 
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We present the /3-expansion of the Helmholtz free energy of the classical XYZ model, with a single- 
ion anisotropy term and in the presence of an external magnetic field, up to order /3 12 . We compare 
our results to the numerical solution of Joyce's [Phys. Rev. Lett. 19, 581 (1967)] expression for the 
thermodynamics of the XXZ classical model, with neither single-ion anisotropy term nor external 
magnetic field. This comparison shows that the derived analytical expansion is valid for intermediate 
temperatures such as kT / J x ~ 0.5. We show that the specific heat and magnetic susceptibility of 
the spin-2 antiferromagnetic chain can be approximated by their respective classical results, up to 
kT I J » 0.8, within an error of 2.5%. In the absence of an external magnetic field, the ferromagnetic 
and antiferromagnetic chains have the same classical Helmholtz free energy. We show how this two 
types of media react to the presence of an external magnetic field. 



1. INTRODUCTION 

The classical limit of quantum spin chains has at- 
tracted attention since this class of models was proposed 
to describe the magnetic interactions of one-dimensional 
models. In 1964 Fishery] calculated the analytical ex- 
pression of the Helmholtz free energy (HFE) of the XXX 
chain in the presence of an external magnetic field. Some 
years later, Joyce included in that model an anisotropy 
term of interaction between the z components of neigh- 
bouring spins, and obtained its exact HFE for a vanishing 
external magnetic field. Since then the thermodynamics 
of the classical XXX chain, in the presence of in-plane 
or transverse magnetic field has been studied at very low 
temperatures (limit of long wavelength) and the mapping 
of this classical spin chain into the Sine-Gordon model|3j 
has permitted one to obtain analytical expressions for 
some termodynamical functions in this region of temper- 
ature. On the other hand, the classical XXZ chain in the 
presence of an external magnetic field in the z direction 
has also been studied numerically 4] . 

Recently we obtained the high temperature expansion 
(HTE) of the HFE of the quantum spin-5 XYZ chainQ, 
with a single-ion anisotropy term and in the presence of 
an external magnetic field, up to order (3 5 . (Here we have 
(3 = ™, where k is the Boltzmann's constant, T is the 
temperature in Kelvin.) The HFE of the classical model 
is obtained from Ref. |5j by taking the limit S — » oo, 
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where S is the spin value (S — 1/2, 1, 3/2, 2, • • • , oo). 
The calculation of this thermodynamical function of the 
quantum model was done by using the cummulant se- 
ries for the one-dimensional models presented in Ref. |(| . 
This approach can also be applied to classical periodic 
chain models with nearest-neighbor interactions. To the 
best of our knowledge, analytical expressions of the HFE 
for intermediate and high temperature regions of the clas- 
sical XYZ model are currently unknown. 

In Ref. 0, we showed that certain thermodynamical 
functions like the magnetization and magnetic suscepti- 
bility per site of the quantum XYZ chain can be well 
approximated, in the intermediate and high temperature 
regions, by their classical results for S > 3/2. The ab- 
sence of results for the quantum XYZ chain for S > 1/2, 
makes the knowledge of the classical model very interest- 
ing in the region of kT > J x , where J x is the coupling 
constant between the closest spins along the x direction. 
Even for other thermodynamical functions, like the spe- 
cific heat per site, the classical models give results with 
the correct order of magnitude for those regions of tem- 
perature. Nowadays the possibility of designing different 
materials with suitable properties turn the analytical ex- 
pressions of the classical thermodynamical functions to 
be very helpful. 

In Ref. [j| we calculated the /3-expansion of the HFE of 
the quantum and classical XYZ model, with single-ion 
anisotropy term and in the presence of an external mag- 
netic field, up to order f3 5 . When applying the method of 
Ref. directly to the classical model we have the advan- 
tage of calculating traces of c-numbers, thus diminishing 
the number of terms to be evaluated. Here is the first 
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time that this method is being applied directly to a clas- 
sical periodic chain with first-neighbour interactions. 

In section 2, we present the Hamiltonian that describes 
the classical XYZ chain with unitary spin (|s*j| = 1) at 
the z-th site. In this model we take into account the pres- 
ence of a single-ion anisotropy term and an external mag- 
netic field along the z direction. In section 3, we present 
the results of Ref.|(| when applied directly to the clas- 
sical XYZ chain; a rule is obtained that allows, in the 
present case, optimization of the algebraic calculations 
of the HFE. In section 4, we compare Joyce's solution^ 
with our results for the particular case of a classical XX Z 
chain {J x — J y ) without the single- ion anisotropy term 
(D = 0) and in the absence of an external magnetic field 
(h — 0), for intermediate and high temperatures. Joyce's 
solution of the classical XX Z chain allows one to obtain 
the behaviour of the specific heat per site at low tem- 
peratures, and then apply Pade's method to extend the 
validity of the HTE to lower temperatures. In section 5 
we study some thermodynamical functions of the classi- 
cal XYZ chain in the presence of an external magnetic 
field. Finally, in section 6 we present a summary of our 
results. 

In appendix^ we present some useful integral results. 
As an example of our HTE of the HFE of the classical 
XYZ chain, we present in appendix iBl its expansion up 
to order /3 6 in the absence of an external magnetic field. 
Its complete expansion, up to order /3 12 , is quite lengthy 
and can be obtained under request to the autors. 



3. THE METHOD OF CUMMULANT SERIES 
APPLIED TO CLASSICAL CHAINS 

In Ref. we presented a survey of the cummulant 
series method and its application to one-dimensional pe- 
riodic quantum chain models with first-neighbour inter- 
actions. Here, we discuss the application of this method 
to any classical one-dimensional chain model subject to 
periodic boundary conditions, spatial translation invari- 
ance, and nearest-neighbor interactions. Following Ref. 
[f|, we obtain the analytical expressions for the HTE of 
the HFE in the thermodynamical limit of such models, 
which can be written as the /3-expansion 
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where trj(lj) equals the dimension of the classical Hilbert 
space of the i-th site, 
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and the auxiliary function tp((3) is given by 
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2. THE CLASSICAL XYZ CHAIN WITH 
UNITARY SPIN 

The classical version of the hamiltonian of the 
anisotropic XYZ chain with unitary spin-5 1 , as shown 
in Eq.(2) of Ref. 0, is 



N 



N 

E 

i=l 



Jy S i S i+l 



{J x sfsf +1 

-hs* + D(s*) 2 }, 



where 



sin(0j) cos( 
cos(0, ; ); 



sin 



(0i) sin(0i), 



(1) 



(2) 



0i and <pi are the polar and azimuthal angles, respec- 
tively, of the classical spin s,; N is the number of sites 
in the periodic chain; h is the external magnetic field 
along the z-axis and D is the single-ion anisotropy pa- 
rameter. The constants J x , J y and J z give the strength 
of first-neighbour interactions between the components 
of the spins. From Q and we define the function H, 
so that 
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The functions H[ nl correspond to the "connected" 
strings with n operators (in the present case, see 

Eq. JJJ) so that m of them are distinct, that is, 
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{M »=1 
The normalized trace is defined as 

(Hj^jj+x . • • Hj m> i m _|_x) = 
tri 1) ... j i m +i(Hj li j 1 _|_i . . . Hi 
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where the indices i\, 



can be equal or distinct. The 



notation tr il: ... iim+ i stands for the trace over the degrees 
of freedom of m + 1 distinct sites in the set: {ii, . . . , i m + 
1}, with m < n. In addition, tXi h represents the trace 
over the degrees of freedom of the ife-th site and li k is its 
identity operator. 

Eq.Q differs from its quantum version on showing a 
normalized trace, rather then a g-trace: in the classical 
case, all H^+i are commuting functions. This property 
greatly simplifies all calculations (see Appendix 0. The 

(3) notation £ " stands for the restriction ^ " Hi = n and 

i=i 
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rii 7^ for i = 1, 2, .., m. The index m satisfies the con- 
dition f ^ to ^ n. Equations @ to (0 are valid for any 
classical one-dimensional chain model subject to periodic 
boundary conditions, spatial translation invariance, and 
nearest-neighbors interactions. In order to reach higher 
orders in (3 in the HTE of the HFE, it is important to 
optimize the calculation of the normalized traces in Eq. 
0. Although the classical model has less terms to be 
calculated than its quantum version, we still have a large 
number of integrals to be done for each order in (3. For 
Hamiltonian it is straightforward to verify that, at 
each site, using the results l|A.4(l and (|A.5|I of appendix 
1X1 that we only have a non-vanishing normalized trace 
((s x ) 1 (s v ) m (s z ) p ) , if and only if I, m and p are all even. 



4. THE /3-EXPANSION OF THE HFE OF THE 
CLASSICAL XYZ CHAIN 

^From the results of Ref. [|| , we can obtain the HTE 
of the HFE of the classical XYZ model, up to order 
/3 5 , by taking the limit of S — » oo in its Eq.(3). By 
using the results of section [21 we obtain here the HTE 
of the HFE, Wdass , up to order (3 12 , for a non- vanishing 
external magnetic field h and in the presence of a single- 
ion anisotropy term in the Hamiltonian • The derived 
expansion is very lengthy and it can be obtained from 
the authors by request. Although this series has a large 
number of terms, it can be easily handled by any CAS 
system (we used Maple). In appendix iBl we present the 
expansion of the HFE of the classical XY Z chain, with 
unitary spin (see Hamiltonian 0, in the absence of an 
external magnetic field (h — 0)), up to order f3 6 . The 
resulting expansion (|B.1|) agrees with the series of this 
thermodynamical function obtained from eq.(3) of Ref. 
0, in the limit S — + oo. 

In this section, we have J x = J y = J and J z — JA, for 
the constants in Hamiltonian for the particular case 
of the classical XX Z model. 



4.1. Comparison of the HTE of the classical XX Z 
chain with Joyce's solution 

Next, we perform some comparisons with known re- 
sults, to verify the interval of temperature where our 
HTEs still hold. Joyce@ obtained the exact solution of 
the classical XX Z chain, without a single-ion anisotropy 
term (D = 0) and in the absence of an external mag- 
netic field (h = 0). In order to check our HTE of the 
classical model and to verify how it extends our previous 
results Q, we compare in Figs, and ^ the Joyce's 
numerical solution for the HFE with our expansions up 
to orders (3 5 and (3 12 . In Figs. ^ and^]d we present 
the respective percental difference between each of those 
expansions and the exact Joyce's solution. In Fig. 
we plot the HFE for A = ±0.5. At first glance of this 
graph suggests that the expansion up to (3 12 is poorer 



than the expansion up to (3 5 . However, a closer look at 
the interval kT/ J = [0.4, 1] (see the detail box) shows 
that the /3 12 -curve coincides with the exact solution up 
to kT/J = 0.53, whereas the /3 5 -curve coincides with it 
only up to kT/J = 0.9. In Fig. [TJa, we take A = ±1 and 
we see that the expansion up to /3 5 coincides with the nu- 
merical solution of Joyce's expression up to kT/J = 0.56 
whereas the series up to f3 12 goes up to kT/J = 0.39, 
which can certainly be considered an intermediate region 
of temperature. 

In Fig. |21 we compare the specific heat per site, 
Cciass(P) = -P 2 -£p? (/3W c ; ass (/3)), obtained from Joyce's 
solution of the classical XXZ chain (with h = and 
D = 0) and the HTE up to /3 6 (see Ref. @) and up to 
[3 13 . We also include (Figs. |2t andJSJi) the respective 
relative percental difference of these expansions to the 
exact solution. In Fig. |2Ji we have J = 1 and A = ±0.5, 
whereas in Fig. |3) we have J = 0.5 and A = ±2. Again, 
in both plots we verify that the /3 13 expansion extends to 
lower temperature the validity of the HTE of the classi- 
cal model in Refs. 0, , as shown in the detail box of 
Fig. [2t. It is simple to derive from Joyce's solution the 
correlation function between the x, y and z components 
of the classical spins of first nearest neighbours. Due to 
the symmetry in the x and y directions in Hamiltonian 
(0, for J x = J y , we have (sfsf+i) = ( s i s i+i)- In Fi S s - 
|3]we plot the correlation function (sf sf +1 ) = dy ^j ass for 
J = 1 and A = 0.5 (Fig. EJi) and for J = 0.5 and A = 2 
(Fig. EJa), derived from the exact result and from our j3- 
expansion (up to order /3 12 ) for the HFE of the classical 
model. From Fig. 03 we see that the /? 12 -curve gives the 
correct maximum of the function (sfsf, x ) for J = 0.5 
and A = 2. 

In Figs. 0|we plot the correlation function (sfsf +1 ) = 
dWciass £ or exaC £ result of the classical model and our 

aJz 

/3-expansion, for the same set of constants J and A as in 
Figs. From Fig. 0k we see that the derived HTE up to 
P 12 gives the correct maximum of this thermodynamical 
function. 

One way to extrapolate the results of our high tem- 
perature series to higher orders in (3 is through the 
Pade approximants (PA), which allows us to combine 
thermodynamical information from both high and low 
temperatures. Among the several approaches to Pade 
approximants!^ , here we employ the two-point Pade 
approximant [lOj to extend our /3-expansion of the spe- 
cific heat per site to low temperatures. 

^From the numerical analysis of Joyce's solution[2j of 
the specific heat per site at very low temperatures, for 
J / 1 and A ^ 1, we realize that it has a polynomial 
behaviour in T. For the range of temperature T £ [0, 0.1], 
it can be chosen to be of fourth order in T, C(T) w 
1 + ayT + a 2 T 2 + a 3 T 3 + a 4 T 4 , where the coefficients a, 
can be adjusted appropriately by linear regression. In 
this case, we have 5 known terms in the region of low 
temperature and 12 known terms in the region of high 
temperature, resulting in a Pade approximant with 17 
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terms. 

FigEJi compares results for the specific heat per site, 
namely the best PAs and Joyce's exact solution, for J = 1 
and A = 0.5. For this constants, we have: a\ = 0.51519, 
a 2 = 1.42213, a 3 = -2.10687 and a 4 = 62.1676. Below 
the graph in[3i, we show the percental deviation of each 
PA with respect to the numerical result of the exact so- 
lution. It turns out that Pn.6 is the best approximation 
of the exact result, within a difference of less than 2% for 
all values of temperature. 

Fig. |SJd shows a similar comparison, but for the pa- 
rameters values J = 0.5 and A = 2. For those constants, 
we obtain numerically: a\ = 0.529029, a 2 — 1.63073, 
a 3 = -8.64101 and a 4 = 86.746. In Fig. Et, we present 
the percental deviation of each PA to the numerical re- 
sult of the exact solution. We see that Pro, 7 is the best 
approximation to the exact result. In the whole interval 
of temperature, its percental difference to the exact solu- 
tion is less than 10% . Although Figs. OJ, and[S|3 refer to 
the same thermodynamical function, each case demands 
specific PAs for the best fitting. 

4.2. The classical model as an approximation to 
the quantum model. 

The classical XYZ chain is a good approxima- 
tion to the quantum chain, for all values of S 
(S = 1/2, 1, 3/2, 2, • • ■ ), in the high temperature region 
(J ZI 8 < 1). In Ref. [g we showed that some thermo- 
dynamical functions of the quantum XXZ chain can be 
well approximated by their classical version for not so 
high temperatures. This region of temperature, where 
the classical and quantum models are equivalent, depends 
on the termodynamical function and on the spin. The 
higher the spin of the quantum model, the more involved 
its numerical solution gets, due to the growing number 
degrees of freedom. However, as far as the XYX model is 
concerned, for S > 2, the quantum HTE of a given ther- 
modynamical quantity can thus be well approximated to 
its classical HTE. 

Some materials are well described by XXZ mod- 
els with higher values of spin. For example, the 
(Ci H 8 N 2 )MnCl 3 , described by a S = 2 model[ll|]; and 
the (CH 3 ) 4 NMnCl3, also known as TMMC, described by 
a5 = 5/2 model[l3ll3t. In Ref. 14], Yamamoto carried 
out Monte Carlo calculations of the thermodynamics of 
the S = 2 XXZ chain with 96 sites, with A = 1, D = 
and h = 0. He obtained the temperature dependence 
of the specific heat per site and the magnetic suscepti- 
bility per site for any temperature. Fig. compares 
the numerical results of Yamamoto Il4j for the specific 
heat per site, for the S = 2 antiferromagnetic case, to 
our /3-expansion of the classical specific heat per site, up 
to order /3 13 . Fig. |BJd shows the relative percental error 
between this two curves, which is less than 2.5% up to 
kT/J = 0.75. In Fig. [7^ we compare Yamamoto's quan- 
tum magnetic susceptibility per site ( X = — ) of the 



5 = 2 antiferromagnetic chain[14j to its equivalent classi- 
cal function obtained from our /3-expansion, up to order 
f3 12 . In the same token, Fig. 0} shows their percental 
difference, which is smaller than 2.5% up to kT/J ~ 0.8. 

In Ref. we verified that the higher the spin, the 
closer quantum and classical thermodynamical functions 
get. The curves, for a given thermodynamical function 
of the quantum models, do not cross for different values 
of spin. Although we do not have a numerical study of 
the thermodynamics of the antiferromagnetic chain with 
S = 5/2, we can affirm that its specific heat and magnetic 
susceptibility per site, for kT/J up to 0.8, differs from 
the classical result in less than 2.5%. This also applies 
to the TMMC [3, since it anisotropy in the z directions 
is small, namely A = 0.016. 



5. THE THERMODYNAMICAL FUNCTIONS 
OF THE CLASSICAL XYZ CHAIN 

^From expression (|B.H of Appendix [0 we verify that 
the HFE of the classical XYZ chain in the absence of 
an external magnetic field (h — 0), is an even function of 
the constants J x , J y and J z , so that this function is the 
same for classical ferromagnetic and antiferromagnetic 
materials with the same single-ion anisotropy D-term. 
As a consequence of this property, we have that some 
thermodynamical functions have well defined parity un- 
der the transformation (J x , J y , J z ) — ► {—J Xl —J y ,—J z )'- 
the specific heat, ((sf) 2 ), the entropy and the mean en- 
ergy are even functions, whereas the first-neighbor corre- 
lation functions, (sfsf +1 ), (sfsf+i) and (s*sf +1 ) are odd 
functions. None of the quantum versions of the above 
mentioned functions have defined parity for h = 5] . 

Figs. |H1 show how each type of medium responds to 
the presence of an external magnetic field. We use the 
fact that the mean energy per site, (£) = j s 
the same for both ferromagnetic and antiferromagnetic 
media at h = 0. Fig. shows the curve of the mean 
energy as a function of J x (3, for J y /J x — 1/3, J z /J x = 
±2/3 and D/J x = -0.6, at h = 0. 

Fig. |Ht> shows, for the ferromagnetic (J z /J x — —2/3) 
and antiferromagnetic (J z /J x — 2/3) cases, the differ- 
ence of the mean energy per site with and without exter- 
nal magnetic field (cf. Fig. 0J>, A£ = (£)\h - (£)\h=o)> 
plotted as a function of h/ J x , for J x j3 = 0.5 and J x j3 = 1. 
Fig. [HJi quantifies this difference, showing the correspond- 
ing percental differences (A£(%) = <£>l f £ 7|f2j, h=0 xl00%). 

The behavior of the mean square of the z-component 
of each spin in the chain at ft = as a function of J x (3, 
(( s f) 2 ) — dw d a ^ BB , can be seen in Fig. |5Jt. The reac- 
tion of the ferromagnetic and antiferromagnetic media 
to the presence of an external magnetic field, ((sf) 2 ) — 
(( s i ) 2 }|/i=0j can be seen in Fig. as a function of h/J x 
at J x (3 = 0.5 and J x (3 — 1. For both graphs, we take the 
same set of values for J y /J x , J z /J x and D/J x as that 
of Fig. [S] for the ferromagnetic and antiferromagnetic 
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media. Fig. ^ is similar to Fig. showing the corre- 
sponding percental difference with respect to ((s| ) 2 }|^ = o- 

In Fig. 1101 we compare the magnetization per site, 
M class = - aW £l a *" , at J x f3 = 0.8, for the ferromag- 
netic case (with J y /J x — 1/3, J z /J x = —2/3, D/J x = 
—0.6) and the antiferromagnetic case (with J y /J x = 1/3, 
Jz/Jx — 2/3, D/J x = —0.6). For these same cases, in 
Fig. ^2 we compare the classic magnetic susceptibility 
per site at h/J x = 0.35. 

Since we are working with unitary classical spins 
(\s,i\ — 1), we define (cos 9i) = (s*, • s^+i), where 9i 
is the angle between the i-th and (i + l)-th spins in 
the chain. In Fig. \12h we plot the function (cos8i) = 

(^-gj- + + ~(fr^j Mclass for the antiferromagnetic case 

(with Jyj J x = 1/3, J z /J x — 2/3, D/J x = —0.6), whereas 
in Fig. 112b we plot (cos 0i) for the ferromagnetic case 
(with J y /J x = 1/3, J z /J x = -2/3, D/J x = -0.6), both 
at h/J x = 0.35. From Fig. U2b we see that at J x (3 = 0.9, 
on the average, neighboring spins are orthogonal to each 
other. For these cases, we also plot in Fig. ^|the corre- 
lation function (sfsf , x ). 

Within the range of the independent variable (h/J x ) 
shown in Fig. 1101 the HTE of the magnetization M is 
almost equal to its exact solution. This range has been 
determined so that the HTEs of leading orders (3 11 and 
j3 12 differs by 0.1% therein. A similar thought guided the 
determination of the ranges of J x (3 in Figs. ^2 an( l H~3l 
regarding the magnetic susceptibility \ an d correlation 
function (sfsf +1 ), respectively. 



6. CONCLUSIONS 

The method developed in Ref. can be equally ap- 
plied to both quantum and classical chains with first- 
neighbor interactions, spatial periodic boundary condi- 
tions and translational invariance. In Ref. p| we cal- 
culated the high temperature expansion (HTE) of the 
Helmholtz free energy (HFE), up to order /3 5 , of the 
quantum spin-S* XYZ chain, with a single-ion anisotropy 
term and in the presence of an external magnetic field. 
From this result, we obtained the HTE of the classical 
version of the model by taking the limit S — > oo. In the 
present paper, we apply the method of Ref. Q directly to 
the classical XYZ (also with a single-ion anisotropy term 
and in the presence of an external magnetic field) thus 
simplifying enormously algebraic calculations. By this 
way, we are able to calculate the /3-expansion of its HFE 
up to order (3 12 . Each coefficient of (3 n in the expansion is 
exact (n = —1, 0, 1, 2, • • ■ , 12). Having a higher order in 
[3 allowed extending the knowledge of the thermodynam- 
ics of the classical XYZ chain up to kT/J x ~ 0.5, which 
might be considered as an intermediate region of temper- 
ature. Joyce's exact result of the classical XX Z chain, 
with no single-ion anisotropy term and no external mag- 
netic field (D = 0, h = 0), permits knowing the behaviour 
of the specific heat per site at very low temperatures. 



This low-temperature information can be combined with 
the high-temperature information by the two-point Pade 
approximants (PA) [ill] . The best PA gives a very good 
description of the classical specific heat per site in the 
whole interval of temperature. 

In the high temperature region (J x f3 <C 1), the quan- 
tum spin-5 XYZ chain (S = 1/2, 1, 3/2, 2, • • •) is well 
described by the classical model. Certainly, the best ap- 
plication of the HFE of the classical XYZ chain is to 
support the study the thermodynamical properties of its 
quantum models for finite spin-5 1 at lower temperatures. 
For S > 2 we have very few numerical analysis of the 
thermodynamical functions of these quantum models due 
to the large number of degrees of freedom to be handled. 
By using the Monte Carlo numerical calculation done by 
Yamamoto|l4j. for a spin-2 antiferromagnetic chain with 
96 sites, we showed that its specific heat and magnetic 
susceptibility per site can be approximated by their re- 
spective classical results, up to kT/ J x w 0.8, within a 
precision of 2.5%. This gives us hope that the thermody- 
namics of the TMCC (S = 5/2) can be approximated by 
the HFE presented in this paper, up to this temperature, 
within a precision higher than 2.5%. 

Finally, from the HTE for the HFE of the classical 
XYZ chain, which can be found in Eq. i|B.l|) of Appendix 
151 we verified that in the absence of an external magnetic 
field the quantum ferromagnetic and antiferromagnetic 
chains have the same classical limit. As a consequence 
of this essentially classical result, some thermodynamical 
functions like: mean energy, ((sf) 2 )}, entropy and the 
correlation functions between first-neighbor spin compo- 
nents have defined parity, at h = 0, under the parame- 
ter transformation [J x , J yi J z ) — > {~J X ,—Jy,—J z )- This 
fact permits studying the reaction of each type of chain 
(ferromagnetic and antiferromagnetic) to the presence of 
an external magnetic field. 
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APPENDIX A: USEFUL INTEGRALS 

The normalized traces of products of operators are ex- 
pressed in terms of surface integrals over unitary spheres, 
each of which represents the state space of a chain site. 
For example, 



(H li2 )= / dhi f d£l 2 H(9i,4>i,02,<h), (A.l) 

J SI JS2 
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where the normalized solid angles d£li are defined 16] as or 



4tt 



■Mi 



(A.2) 



and the function TL has been defined in Eq.J3J. We have 
9i G [0, 7r] and cpi G [0, 2tt]. Another example is 



(H 12 H 2 ,3) 



dni / d« 2 / <m 3 ■H{e ll (j) 1 ,e 2l (p 2 )- > 



SI ^ S2 J S3 

x n(e 2 ,<h,03,<h)- 



(A.3) 



The generalization of these formulas is straightforward. 
Because of the structure of the Hamiltonian and the 
expression of spin operators 10, the integrals related to 
the i-th site which contribute to -f/j 1 ^ (see Eq. Q), turn 
out to be either 



dip (sin (/>) m (cos < 



1 TV m+1 'ipf 't+l ^ 

^(i+(-in(i +(-!)-) y ( Un^ ■ 



(A.5) 



where n, m = 0, 1, 2, . . . and T is the gamma function. 



/I m ' n) = / d9 (sin0) m (cos0) r 



i(i+(-in(i +(-in y r iy +2 ] } 

+ 2 l 1 J J n+1 r(2ttf±2) ' 



(A.4) 



APPENDIX B: THE /3-EXPANSION OF THE 
HELMHOLTZ FREE ENERGY IN THE ABSENCE 
OF A MAGNETIC FIELD 



The high temperature expansion of the HFE of the 
classical XYZ (S — > oo) in the absence of an external 
magnetic field (Hamiltonian with h = 0), up to order 
P 6 , is 
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(B.l) 



In the absence of an external magnetic field (h — 0), and J z of the Hamiltonian Q. We have confirmed this 
this is an even function of the coupling constants J x , J y property up to order /3 12 . 
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FIG. 1: (a) For A = ±0.5, comparison of the numerical solution of Joyce's expression (solid line) and our HTEs up to order 
/3 5 (dashed line) and order /3 12 (dashed-dotted line), for the HFE of the classical XXZ chain (Wdass), as a function of kT/J. 
In (c), the relative percental differences (AW(%)), with respect to Joyce's exact solution, of the /3 5 (dashed line) and the (3 12 
(dashed-dotted line) expansions, also as functions of kT/J. In (b) and (d) a similar comparison is performed, for A = ±1. 
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FIG. 2: (a) Specific heat per site Cdass of the classical XX Z chain, as a function of kT, T being the absolute temperature, for 
J = 1 and A = ±0.5. The upper pane shows Joyce's solution (solid line), and our expansions up to /3 6 (dotted-dashed line)|f| 
and up to f3 ls (dashed line). The lower pane (c) shows the percental differences (AC(%)), with respect to Joyce's solution, 
of the /3 6 (dotted-dashed line) and /3 13 (dashed line) expansions, also as functions of kT. Panes (b) and (d) show a similar 
comparison, for J = 0.5 and A = ±2. 
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FIG. 3: The correlation function (sf sf +1 ) between the classical spin a;- components of first neighbouring sites. The solid line 
stands for the exact function, whereas the dashed line stands for the HTE up to order /3 12 . In (a) we have J = 1 and A = 0.5; 
in (b) we have J = 0.5 and A = 2. 




FIG. 4: The correlation function (sfsf +1 ) between the classical spin ^-components of first neighbouring sites. The solid line 
stands for the exact function, whereas the dashed line stands for the HTE up to order (3 12 . In (a) we have J = 1 and A = 0.5; 
in (b) we have J = 0.5 and A = 2. 
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FIG. 5: Upper panes (a) and (b) compare exact classical specific heat per site and their best PAs. Lower panes (c) and (d) 
show the percental difference of each PA with respect to the exact result. In (a) and (c) we have J = 1 and A = 0.5, whereas 
in (b) and (d) we have J = 0.5 and A = 2. 




FIG. 6: (a) The dots correspond to the Monte Carlo (MC) calculation of the specific heat per site for the S = 2 antiferromagnetic 
chain jJJ, for D = and h = 0, as a function of kT, where T is the absolute temperature. The solid line corresponds to our 
HTE of the classical specific heat up to order /3 13 . (b) The relative percental difference of the two curves. 
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kT/J 



FIG. 7: In (a) we compare Yamamoto's MC calculation of the magnetic susceptibility per site of the S = 2 antiferromagnetic 
chain |l4| (dots and dashed line) and the HTE of its classical equivalent, up to order f3 12 (solid line). In (b), we have the 
percental difference of the two functions plotted in (a). We have A = 1, D = and h = 0. 
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FIG. 8: (a) The mean energy per site (e) of ferromagnetic (J z /J x = —2/3) and antifcrromagnetic (J z /J x = 2/3) media in the 
absence of an external magnetic field, (b) The difference Ae of the mean energy per site at finite h and its value at h — versus 
h/J x . We plot the antiferromagnetic case for J x f3 — 0.5 (solid line) and J x /3 = 1 (dotted line); and also the ferromagnetic 
case for J x f3 = 0.5 (dashed line) and J x f3 = 1 (dotted-dashed line). In (c) we have the correspondent percental differences 
Ae(%)with respect to the h = case; the same graphical convention for lines is used. In all figures we have J y /J x = 1/3 and 
D/J x = -0.6. 
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FIG. 9: (a) The mean value of the squared z-component of spin per site ((sf) 2 ) for ferromagnetic (J z /J x = —2/3) and 
antiferromagnetic (J z /J x = 2/3) media is the same, in the absence of an external magnetic field, (b) The difference of ((sf) 2 ) 
at finite h and its value at h — versus h/J x . We plot the antiferromagnetic case at J x f3 — 0.5 (solid line) and J x f3 — 1 
(dotted line); and the ferromagnetic case at J x f3 = 0.5 (dashed line) and J x /3 = 1 (dotted-dashed line). Fig. (c) shows the 
corresponding percental differences with respect to the h = case; the same graphical convention for lines is used. In all figures 
we have J y /J x = 1/3 and D/J x = —0.6. 
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FIG. 11: The classical magnetic susceptibility per site Xciass as a function of J x f3 with J y /J x = 1/3 and D/J x = —0.6 at 
h/J x = 0.35, for the ferromagnetic (J z /J x = —2/3, dashed line) and antiferromagnetic (J z /J x = 2/3, solid line) cases. 
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J x (3 J x (3 

FIG. 12: The function (cos ft) as a function of J x f3, for i = 1, 2, • • • , N with J y /J x = 1/3 and D/J x = -0.6 at h/J x = 0.35. 
Fig. (a) shows the antiferromagnetic case (J z / J x = 2/3) and Fig. (b) the ferromagnetic case (J z /J x = —2/3). 




FIG. 13: The correlation function (sfsf +1 ) between the spin z-components of first neighbors as a function of J x f3. The dashed 
and solid lines describe the ferromagnetic case (J z /J x = —2/3) and antiferromagnetic case (J z /J x = 2/3), respectively. We 
take J y /J x = 1/3 and D/J x — —0.6 at h/J x = 0.35. 



